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Constraining the star formation rate in the Solar 
neighbourhood with star clusters 
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ABSTRACT 

This paper investigates the star formation rate (SFR) in the Solar neighbourhood. 
First, we build the local age distribution function (ADF) with an updated sample of 
442 star clusters located at less than 1 kpc from the Sun. Next, we define the SFR, 
compute the individual mass evolution of a population of artificial clusters covering 
the broad range of parameters observed in actual clusters, and assume 100 M Q as the 
low-mass limit for effective cluster observation. This leads to a simulated ADF, which 
is compared to the low-noise Solar neighbourhood ADF. The best match corresponds 
to a non-constant SFR presenting two conspicuous excesses for ages < 9 Myr a nd 
between 220 — 600 Myr (the local starburst). The average formation rate is SFR w 
(2500 ± 500) M Q Myr" 1 , corresponding to the average surface formation rate T^sfr ~ 
(790 ± 160) M Myr -1 kpc -2 . These values are consistent with the formation rate 
inferred from embedded clusters (ECs), but much lower (< 16%) than that implied 
by field sta rs. Bo th the local starburst and the recent star formation period require 
SFR ~ 2 x SFR to be described. The simulations show that 9 1.2 ±2. 7% of the clusters 
created in the Solar neighbourhood do not survive the first 10 Myr, which is consistent 
with the rate of EC dissolution. 

Key words: ( Galaxy:) open clusters and associations: general 



INTRODUCTION 



It is currently accepted that star formation is scale-free and 
hierarchical, with high-velocity turbulent gas forming large- 
scale str uctures, while lo w-velocity compression forms small 
clumps (Elmegreen 2008). A natural consequence of this sce- 
nario is that young stellar groupings would be hierarchically 
clustered, with the great star complexes at the largest scales 
and the OB associations and subg roups, clusters and cluster 
sub-clumps, at the smallest (e.g. lEfremovl ll995). The bot- 
tom line is that star formation appears to occur primarily 
in clusters. Consequently, star clusters appear to be excel- 
lent tracers of the star formation rate (SFR) and history of 
the Galaxy, provided that their fundamental parameters - 
especially the age - are known. 

However, star clusters have limited lifetimes. At the ear- 
liest phases, most of the embedded star clusters (ECs) dis- 
solve into the field on a time-scale of a few 10 7 Myr. Disso- 
lution here occurs mainly because the gravitational poten- 
tial can be rapidly reduced by the impulsive gas removal 
by supernovae and massive star winds associated with this 
early period. Thus, a significant fraction of the stars, espe- 
cially of low mass, end up moving faster than the scaled- 



down escape velocity, an d may escape to the field (e.g. 
I Goodwin fc Bastianll200o1 ). 

Dissolution also affects the open clusters (OCs) that 
survive the early phase. Even around the Solar circle, most 
OCs dissolve long b efore reaching an age of ~ 1 Gyr (e.g. 
iBonatto et alj |2006). The disruption time scale, in general, 
scale s with mass as £<jm ~ M 0S2 (e.g. lLamers fc Gielea 
2006), which means that Solar neighbourhood clusters with 
mass in the range 10 2 — 10 3 Mq dissolve on a timescale 
75 < tdis(Myr) < 300. This occurs because OCs contin- 
ually undergo mass segregation and evaporation, tidal in- 
teractions with Galactic substructures, shocks with giant 
molecular clouds (GMCs), as well as mass loss due to stel- 
lar evolution. By decreasing the total cluster mass - and the 
collective gravitational potential, these processes affect the 
internal dynamics and accelerate the cluster dynamical evo- 
lution. Eventually, the majority of the OCs dissolve in the 
Galactic stellar field. 

Thus, the only way for recovering the SFR using the 
observed star cluster population is by taking the mass- 
loss processes affecting clusters of different masses and 
orbits into account. To this effect, semi-analytical de- 
scriptions of the ma in mass-loss processes have becom e 
available recently (e.g. lLamers. Baumgardt fc Gielesl (|201Ct ): 



2 C. Bonatto and E. Bica 



iKhalisi, Amaro-Seoane fc Spurzemll2007l h In the present pa- 
per we employ these process to simulate the mass evolution 
of a population of artificial clusters. Properties of the local 
SFR are then derived by comparing the simulated age dis- 
tribution function (ADF) with that built for clusters in the 
Solar neighbourhood. 

This paper is organised as follows. In Sect. [2] we build 
the ADF for the Solar neighbourhood. In Sect. [3] we briefly 
discuss the cluster mass-loss process. In Sect. 0] we simulate 
the observed ADF and use it to constrain the local SFR. 
Concluding remarks are given in Sect. [5] 



2 BUILDING THE OBSERVED CLUSTER ADF 

The number of star clusters with accurate age and dis- 
tance determinations has been steadily increasing over 
the last years. This is particularly true for the very 
young clusters, most of which so embedded in their 
parent gas and dust cloud that their stellar content is 
essentially inacessible to optical photometry. However, 
the availability of uniform, wide-field, and rather deep 
near-infrared surveys (e.g. 2MASS3, GLIMPSE0, some 
covering essentially all the sky (2MASS), has led to the 
discovery - and allowed a robust parameter derivation - 
of m a ny such embedde d clusters (ECs - e.g. iDutra et aF 



20031; iBica et al.1 |200£ 


; iKumar. Kamath & Davis| |2004); 


Kumar, Keto & Clerkii 


12006); Bonatto. Santos Jr. & Bical 


20061; Ortolani et al.l 


20081; Bonatto & Bical 2009b; 


Bonatto & Bical l2009d; iBonatto & Bical l2010al). to- 



gether with some old and /or very redd ened PCs (e.g 



Froebrich. Scholz fc Raftervl 12007 



20071; IBonatto fc Bical l2008t 



.Bonatto fc Bica 
Bonatto fc Bical l2009ah 



Froebrich et alj|20ld : IBonatto fc Bicall2010bl )" 

Together with the several hundred objects already in- 
dexed in the widely-used star cluster databases WEBDA0 
and DAMLO50, the recent discoveries (together with param- 
eter derivation for poorly-studied and/or unstudied objects) 
are proving invaluable in constructing a more detailed pic- 
ture of the cluster ADF, especially in the Solar neighbour- 
hood and for very young clusters (Sect.|4)|. 

We started by searching WEBDA and DAML02 for 
clusters with available age (4a) and distance from the Sun 
(cIq). However, given the amount of new data routinely 
published, it usually takes a considerable time for both 
databases to incorporate the recently discovered clusters, 
or to update previous entries with newly derived parame- 
ters. Then, we complemented the sample by searching the 
recent literature for clusters that still are not listed in ei- 
ther database. Cluster designations and coordinates have 
been checked among all sources (WEBDA, DAML02, and 
literature) to avoid duplicity. When multiple values of age 
and/or distance occurred, we adopted those based on colour- 
magnitude diagrams (CMDs), or the more recent. The fi- 



1 The Two Micron All Sky Survey, All Sky data release 
llSkrutskie et al.ll2006h 

2 Galactic Legacy In frared Mid-Plane Survey Extraordinaire 
llBeniamin et al.ll2003f) 

3 www.univie.ac.at/webda 

4 Catalog of Op tically Visible Open Clusters and Candidates, 



nal sample contains 1718 clusters (ECs and OCs) with age 
and distance, of which 442 are closer than 1 kpc from the 
Surjf]. By far, most of the parameters have been taken from 
WEBDA and DAML02, which do not provide measurement 
uncertainties. Thus, based on our experience in working 
with clusters of different ages and distances, we adopted the 
following uniform error attribution: 10% for d© < 1 kpc, 
15% for 1 < d©(kpc) < 5, 20% for 5 < d©(kpc) < 9, 
and 25% for d© > 9 kpc; 35% for t A < 20Myr, 30% for 
20 < U(Myr) < 100, 20% for 100 < i^(Myr) < 2000, and 
50% for t A > 2000 Myr. 

Age uncertainties are explicitly incorporated into the 
ADF, which is defined as the fractional number of clus- 
ters per Myr, ADF = dN/dtA- Formally, if measurements 
of a given parameter \ are normally (i.e. Gaussian) dis- 
tributed around the average x with a standard deviation a, 
the probability of finding it at a specific value x i s given by 

P( x ) = -/=- e"^ 2 ^) 2 . To build this ADF we first define 
a set of bins spanning the whole range of ages, and having 
widths that increase with age (to account for the decreas- 
ing number of clusters at older ages). Then, for a cluster 
with age and uncertainty tx±ff, we compute the probabil- 
ity that the age corresponds to a given bin, which is simply 
the difference of the error functions at the bin borders. By 
doing this for all clusters and age bins, we have the number- 
density of clusters in each age bin. By definition, the integral 
of the ADF over the whole range of ages is the number of 
clusters. Subsequently, bin widths can be adjusted to min- 
imise the errors, so that the resulting ADF has statistically 
meaningful values over all ages. 



3 OVERVIEW OF THE MASS-LOSS 
PROCESSES 

Star clusters lose mass continually by a combination of pro- 
cesses associated with stellar evolution and dynamical inter- 
actions (both internal and external to the cluster). Robust 
analytical descriptions of the mass-loss processes - for clus- 
ters characterised by a wide variety of parameters and or- 
biting in different environments - have become available in 
recent years, wi th model parameters derived from theoreti- 
cal gr ounds fe.g. lSpit zcr 1987; Lamers. Baumgardt fc Gielesl 

ma 



2010h and N-body simulations (e.g. iBaumgardt fc Makinc 
2003; iGieles fc Baumgardtll2008h . Formally, the time-rate of 
change of mass of a cluster that was formed with the mass 
Mi = M (0) can be expressed as: 



dM 
~df 



E 

p=i 



(dM_\ 
\~dTJ, 



(1) 



IDias et al.l 



ot Op 
j2002h 



where the mass- loss process p are: (1) stellar evolution, 
(2) tidal effects by a steady field, (3) shocks with spi- 
ral arms, (4) encounters with GMCs, (5) evaporation, 
and (6) ejection. For proce sses (1) - (4) w e adop t the 
semi-analytical approach of ILamers fc Gielesl (|2006t ) and 
ILamers. Baumgardt fc Gielesl ( 2010T ). In what follows we as- 
sume that masses are always expressed in Solar masses (M©) 



5 Upon request, the table with the age and distance from the Sun 
for the 442 nearby clusters is available from one of us (C.B.). 
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and time in Myr; also, we write the equations in terms of the 
remaining- mass fraction jj, = M/Mi. For clarity, we provide 
below a brief description of these processes. 



3.1 Stellar evolution 

Since we are dealing essentially with disk clus- 
ters, we adopt the Solar-me ta Uicity approxima- 
tion given by lLamers fc Gielesl (|2006l ). which is 
based on single-stellar popu l ation GALEV models 
jAnders fc Fritze-v. Alvenslebenl [2003) and a Salpeter- 



like (|Salpetej|l955h mass function: 



(f) --(IK- 1 "" 



s(ii)+qq a +b] 



with qq = log(t) - 1.0, a = 0.255, and b = -1.805. This 
process is relevant for t > 10 Myr. 



3.2 Tidal effects 

We adopt the semi-analytical description of the 

mass-loss of clusters on a steady-tidal field of 

lLamers. Baumgardt fc Gielesl l|2010f ). which is expected 
to apply to a broad range of cluster conditions and orbit 
environments: 



\dt J tidal 



toW 



with t — 



tR (kz*) ( R o ) ( 



22okm s 



t R = 13.3 Myr 



for clusters with an initial concentration factor of the density 
King profile W = 5, and t R = 3.5 Myr for Wo = 7; 7 = 0.65 
for Wo = 5 and 7 = 0.8 for Wo = 7, thus, both 7 and tn 
are not independent parameters, since they can be expressed 
as a function of Wo; £ is the cluster orbit eccentricity; the 
average stellar mass at time t is given by the interpolation 
formula m(t) = 0.6193 - 0.0362r - 0.01481r 2 + 0.0022-r 3 , 
where r = log(t); Ra is the cluster's Galactocentric dis- 
tance, and vg is the (assumed constant) rotation velocity. 
In what follows we take vg = 200 km s" 1 for orbits around 
the Solar neighbourhood. 



3.3 Shocks with spiral arms 

The energy gain and the mass loss due to the disruptive 
effect of spiral-arm passages of star clusters with planar and 
circular orbits arou nd the centres of galaxies has been thor- 
oughly studied by iGieles. Athanassoula fc Portegies ZwartJ 
J200 



I 



In particular, for the 



Gieles. Athanassoula fc Portegies ZwartJ 
Lamers fc Gielesl (|2006h find: 



Solar n eighb ourhood, 
(|2007l l and 



\dtJ, 



-5 x 10" 



(Mi/1O 4 M ) - 7 ' 



3.4 Encounters with giant molecular clouds 

Similarly to the spiral arms, lLamers fc Gielesl i|2006h esti- 
mate the energy gain and mass loss due to encounters be- 



tween clusters and GMCs. Assuming GMC parameters typ- 
ical of the Solar neighbourhood, they find: 



\dt). 



= -5 x 10" 



(M l /10 4 M s )0-7- 



The mass loss due to GMCs is thus 10 times higher than 
that of the spiral arms. 



3.5 Ejection and evaporation 

Stars can also escape from a cluster by means of ejection 
and evaporation. Ejection occurs when, after a single, close 
encounter with another star, a member ends up with excess 
velocity with respect to the escape velocity. Evaporation, 
on the other hand, is related to a series of more distant 
encounters that gradually increase a star's energy, eventually 
leading it to escape from the cluster. With the respective 
time scales taken from Khalisi, Amaro-Seoane & Spurzem 
(2007, and references therein), the mass loss associated with 
both processes are, respectively: 



and 



\ dt J el 



m(Gfi) 1/2 
ej - 46(M^3 m )V2' 

fh(Gn) 1/2 ln(7 c A0 



13.8(M^| m )V2 ' 

where the half-mass radius is given by Rhm = 
3.75 (M/10 4 ) 01 l|Larsen|[200i ). -y c = 0.11 is the Coulomb 
factor, N = M(t)/fh(t) is the number of stars at time t, and 
G is the gravitational constant. 

3.6 The adopted procedure 

In summary, given a cluster of initial mass Mi, and a set 
of model parameters (Rgc, Wo, 7, e, t R ), Eq. [T]should be 
solved to find the mass still bound to the cluster at a later 
time t, M = M(t). 



4 INTERPRETING THE LOCAL ADF 

The Solar neighbourhood ADF has been subje c t to p revi- 
ous investigation . In p articular, lLamers et all (|2005h and 
lLamers fc Gielesl (|2006l ) employed the cluster dissolution 
processes 1—4 (Sect. |3} to study the ADF of 114 PCs 
closer than <i© = 0.6 kpc taken from iKharchenko et ail 
(2005). Their ADF can be reasonably well described by 
a nearly-constant SFR (in bound clusters of mass 10 2 < 
M(M ) < 3 x 10 4 ) of SFR = 400 M Myr" 1 , corresponding 
to a surface formation rate T,sfr ~ 350 M Myr" 1 kpc -2 , 
which is considerably lower than the Y>s fr = 700 — 
1000 M Myr" 1 kpc" 2 inferred from ECs (|Lada fc Ladal 
l2003h . and the S SFfl = 3000 - 7 000 M Myr" 1 kpc" 2 from 
field stars (| Miller fc Scald Il979l l. They also find evidence 
of a local burst of star formation that took place betw een 
250— 600 Myr ago (also present in lPiskunov et al.ll2006h . 

Besides cluster dissolution, observational completeness 
is also important for the ADF shape. High absorption and 
crowding, usually associated with fields dominated by disc 
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and bulge stars, affect cluster detec tability, especia l ly the 
faint and/or poorly-populated ones. iBonatto et al.l (|2006l ) 
show that most of the intrinsically faint and/or distant 
clusters must be undetected in the field, particularly in 
bulge/disc directions. They also find that the completeness- 
safe zone might extend up to d « 1.4 kpc inside the Solar 
circle, but considerably more outside. 

Thus, to minimise observational completeness effects, 
we restrict the analysis to the region d,Q ^ 1 kpc, which 
reduces the number of clusters to 442, but still a statisti- 
cally significant sample. Besides, this condition also guar- 
antees that we are dealing with clusters subject to some- 
what similar physical conditions. By construction (Sect. [2]), 
the resulting ADF (Fig. [1} is rather smooth, with small er- 
rors over all ages, except for tA > 3 Gyr, which reflects the 
scarcity of Gyr-old clusters in the region. Also conspicuous 
are the excesses between 200 < ^(Myr) < 600, and es- 
pecially for tA < lOMyr (see bel ow). Since we work with 
a significantly larger sa mple than lLamers et all (|2005l ) and 
lLamers fc Gielesl (|2006l ) , our approach provides a better def- 
inition of structures in the ADF. 

4.1 ADF simulation 

Instead of assuming particular values for the several model 
parameters (Rgc, Wo, 7, ta, e, number, age, initial mass, 
and the orbit of clusters), we adopt a simplifying approach. 
Starting with a pre-defined SFR, we compute the time- 
evolution of the individual mass of star clusters created with 
parameters cov ering the observed range in orbits around 
the Solar circle (|Lamers. Baumgardt fc GieledlioibT ). We as- 
sume that the relevant model parameters can take on any 
value within 0.0 ^ e ^ 0.8 (also allowing for circular orbits), 

5 s: Wo ^ 7, 0.65 sC 7 0.8, and 3.5 t R (Myr) s; 13.3 
(the latter two are expressed as a function of Wo as in = 
37.8 - 4.9H/ and 7 = 0.275 + 0.075W ); also, given the re- 
striction rig ^ 1 kpc, we have 7.5 ^ i?Gc(kpc) ^ 9.5. Each 
cluster is created (at time tA) with an individual mass in the 
range Mmin < Mi(Mg) < M max , but following t he number- 
frequ ency (or probability) dN/dMi oc M~ 2 (e.g. lElmegreenl 
;2OOi£0. We adopt M mi „ = 10 M Q and M max = 7.5 x 10 4 M 
as the minimum and maximum initial cluster mass, respec- 
tively. Then we solve Eq. [T] to find the mass remaining in 
each cluster after evolving for a time fro m t = to t = tA, 
and b uild the simulated ADF. Similarly to lLamers" fc Gielesl 
(2006), we assume 100 Mq as the minimum cluster mass for 
completeness not being an issue to detectability, i.e, only 
clusters more massive than 100 Mq are accounted for in the 
ADF. According to this approach, the number of free pa- 
rameters reduces essentially to the SFR shape. 

The following steps are taken: (i) We assume that all 
star formation occurs in clusters, define the SFR and com- 
pute the total mass [Mtot = / SFR(t) dtj converted into 
stars over the Galaxy age; in the case of a non-constant 
SFR, this should be done separately for each SFR segment 
(or period of time), (ii) We assign the initial mass Mi (see 
above) for one cluster. (Hi) Next, we randomly assign the 
cluster age tA (it should be within the corresponding SFR 

6 This is done by Mi = i_ n (i_M m *™ /m )' wnere n is uni- 
formly distributed within [0,1]. 



time segment), and the model parameters, (iv) Steps (ii) and 
(in) are repeated for more clusters, until the sum of initial 
masses equals M to t- (v) Then we solve Eq. [T]for all clusters 
from t — to tA, and keep only those with M{tA) > 100 Mq. 
(vi) Finally, stochastic effects are minimised by repeating 
steps (ii) to (v) several times (N S i m ). In the present case we 
use N a i m — 100, and solve Eq. [T]for M = M(t) by means of 
a 4 th order Runge-Kutta method, with a fixed time step of 
lMyr. 

4.2 Constraining the local SFR 

Given the conspicuous excesses (Fig. [TJ, it is clear that 
the Solar neighbourhood cluster ADF cannot result from 
a constant SFR. Instead, the best match has been ob- 
tained with the segmented SFR: SFR = 5040 M Q Myr -1 
(Zsfr ~ 1600 MoMyr" 1 kpc" 2 ) for t A «S 9Myr, SFR = 
48OOM Myr- 1 (Zsfr ~ 1500 M Q Myr -1 kpc" 2 ) for 220 < 
t A (Myr) < 600, and SFR = 2400 M Q Myr" 1 (E SFJ ? w 
760 M Myr" 1 kpc -2 ) elsewhere (bottom panel of Fig. [Q. 
Integrating the segmented-SFR ADF over time yields the 
same number of clusters as that of the observed ADF. Also, 
the lcr boundaries of the adopted solution (shaded region 
in Fig. [T]), obtained after 100 different simulations, are very 
narrow over most of the cluster ages. To check how con- 
strained is this solution, we apply the same procedure to 
SFRs that are 20% higher and lower than the adopted SFR 
over all time periods. Overall, both SFRs provide a much 
poorer description of the observed ADF than the segmented 
SFR, especially for clusters younger than ~ 300 Myr. Thus, 
we adopt 20% as the SFR uncertainty. As expected, the 
ADF produced by a constant SFR (2500 Mq Myr" 1 ) fails to 
describe the observed excesses. 

Considering the above arguments, the average values 
(over the full age range) of the SFR and surface SFR 
are S~FR w (2500 ± 500) M Q Myr" 1 and Esfa » (790 ± 
160) Mq Myr -1 kpc -2 . Given the higher number of clus- 
te rs in our samp l e, we find a T*sfr more than twic e that 
of lLamers et all |2005l) and lLamers fc Gielesl |2006l ) , con- 
siderably more cons istent with that inferred from ECs by 
lLada fc Ladal (|2003l ). but still signifi cantly lower (~ 16% ) 
than the rate implied by field stars l|Miller fc Scald fl979l ). 
We also compute the fraction of clusters that do not sur- 
vive the first 10 Myr, 91.2 ± 2.7%, which agrees with the 
rate of cluster diss olution after the embedded phase (e.g. 
lLada fc Ladall2003l ). 

On average, ~ 3.3 x 10 5 clusters are created (at dg < 
1 kpc) in each simulation with the segmented SFR. This 
implies that less than 0.15% of the clusters ever formed 
in the Solar neighbourhood can be observed (i.e, having 
mass > 100 Mq). And, if the same SFR applies to Rgc = 
7.5 — 9.5 kpc, ~ 10 7 clusters have been created in this region 
since Galaxy formation. 



5 SUMMARY AND CONCLUSIONS 

In this paper we analyse properties of the Solar neighbour- 
hood SFR, by means of a statistically significant ADF built 
with 442 star (embedded and open) clusters closer than 
1 kpc from the Sun. By adopting a simplifying approach, 
in which the mass evolution of artificial clusters is followed 




Figure 1. Top: The Solar neighbourhood ADF (circles) is well reproduced by a segmented SFR (heavy-solid line) with enhanced rates at 
^ 9Myr and 220 — 600 Myr and small Stochastic fluctuations over 100 simulations (shaded region). Also shown are the ADFs produced 
by Hsfji's 20% higher (dashed line) and lower (dotted) than above, together with a constant Y,gpR (heavy dot-dashed line). Inset: 
Blow-up of the local starburst period. Bottom: The Hsfr's (in Mq Myr -1 kpc -2 ) used in the top panel. 
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over time, we reduce the problem to essentially finding the 
SFR. 

The artificial clusters embody parameters and condi- 
tions expected to apply to most actual star clusters orbiting 
not far from the Solar circle. To simulate the observed ADF, 
we employ semi-analytical descriptions of the mass-loss pro- 
cess responsible for cluster dissolution, and assume that only 
clusters with a present-day mass above 100 M can be ef- 
fectively observed (i.e., take part in the ADF). 

The best match between observed and simulated ADFs 
corresponds to a non-constant SFR, with enhanced rates for 
ages the < 9 Myr and 220 — 600 Myr (the so-called local star- 
burst). The average rate is ~SFR « (2500 ± 500) M Myr -1 , 
corresponding to the average density T,sfr ~ (790 ± 
16O)M0Myr _1 kpc" 2 . These values agree with the forma- 
tion rate inferred from ECs, but represent only ~ 16% of the 
rate implied by field stars. Both the local starburst and the 
recent formation (^ 9 Myr) have ADF amplitudes suggest- 
ing periods with a SFR about twice the average value. We 
also find that 91.2 ±2.7% of the clusters created in the Solar 
neighbourhood dissolve before 10 Myr, which is consistent 
with the rate of EC dissolution. 
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